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ABSTRACT 
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The goals of the IRIS Experiment are the measurement of 
parameters describing the vertical structure of the atmosphere 
such as temperature, humidity and cloud height. This information 
is extracted from the thermal emission spectrum between 5 and 
20 At. A Michelson type interferometer capable of resolving to 
5 cm -1 is described in detail. //). , 
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THE INFRARED INTERFEROMETER 
SPECTROMETER EXPERIMENT (IRIS) 

(VOLUME H - METEOROLOGICAL MISSION) 


SUMMARY 

Meteorological satellite experiments have mapped the earth in several wave- 
length regions of the electromagnetic spectrum. Two dimensional maps were 
generated which proved very useful for research as well as operational pur- 
poses. However, the most promising tool of modern meteorology, numerical 
forecasting by high-speed computers, must treat the atmosphere in three di- 
mensions. Atmospheric processes in the troposphere are governed primarily 
by hydrodynamic forces. Integration of the Navier-Stokes equation or simplified 
forms of it require wind velocities, temperature profiles, water vapor distribu- 
tions, cloudiness and surface pressure. Radiative heating and cooling are less 
important for short term forecasting, but influence the atmospheric processes 
over longer periods strongly. Today, only the surface temperature and cloudi- 
ness are under satellite surveillance. 

The purpose of the IRIS experiment is to provide the temperature profile 
and the water vapor distribution. Wind velocities and surface pressure are still 
beyond the capabilities of remote measurements, but the potential to measure 
the surface pressure exists in measurements of the thermal emission spectrum. 
Besides the aforementioned goals, various other tasks can be accomplished by 
the IRIS experiment: ozone distributions, surface temperatures, and cloud height 
estimates may be derived from the data. 

The above information will be extracted from the thermal emission spec- 
trum within the interval of 2000 cm" 1 (5 microns) to 500 cm" 1 (20 microns) and 
possibly to 300 cm" 1 (33 microns) with a spectral resolution of 5 cm -1 . The 
IRIS experiment is deliberately aimed at a wide spectral range to permit cover- 
age of absorption regions of water vapor, carbon dioxide and ozone as well as 
the atmospheric windows. 

Design considerations for the IRIS experiment are analyzed in detail. Spec- 
tral resolution, spectral range, and spatial resolution were influenced strongly 
by requirements to make the results available to numerical forecasting. 
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For Nimbus, the 8-degree field of view of the instrument covers a circle 
of about 160 km on the earth's surface, with the optical axis pointed toward the 
subsatellite point. Each spectrum is taken within 11 seconds. Image motion 
compensation is provided to eliminate smear. Calibration in orbit assures ac- 
curacy of measurement. 

A Michelson-type interferometer is used in this experiment. A breadboard 
model of the instrument exists and is presently undergoing calibration and en- 
vironmental tests. This report includes the design criteria for this interfer- 
ometer and a detailed analysis of the achievable signal-to-noise properties. 

The data reduction of the IRIS experiments requires a power spectrum 
analysis. Computer programs for this task were constructed and samples of 
simulated interferograms and analyzed spectra are shown. 

Besides the application of the IRIS experiment to meteorology on earth, the 
same experiment and essentially the same instrumentation are capable of pro- 
viding unique information on the atmospheres of Venus, Mars, Jupiter. 


ix 
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THE INFRARED INTERFEROMETER 
SPECTROMETER EXPERIMENT (IRIS) 

(VOLUME H - METEOROLOGICAL MISSION) 


INTRODUCTION 

Scientific experiments have been performed from meteorological satellites 
since 1959. Vanguard II, Explorer VII, the TIROS series, and finally Nimbus 
are milestones on this road of increasing success. Mapping of the earth by tele- 
vision cameras and by infrared radiometers allows recognition of cloud patterns, 
identification of frontal systems, and tracking of major storms. The measure- 
ment of reflected solar radiation and thermal radiation emitted by the planet as 
a function of latitude and season provides data for the study of the heat balance 
and of climatology. Comparison of equivalent blackbody temperatures of the 
clouds as measured in the atmospheric windows (8-12 p. or 3. 5-4.1 p) with clima- 
tological temperature profiles makes it possible to estimate the height of the 
upper limits of clouds. These early estimates and the more recent measure- 
ments of TIROS 7 in the I5p C0 2 band were the first attempts to derive atmo- 
spheric parameters (cloudiness and temperatures, in this case) as a function of 
height. In general, however, radiometric mapping with relatively high spatial, 
but low spectral, resolution yields primarily cloud or moisture patterns in a two- 
dimensional sense but provides very little information on the vertical distribu- 
tion of temperature, water vapor, or other parameters. Atmospheric motions, 
the vertical transport of moist air, and other meteorologically interesting phe- 
nomena can be recognized, described, and finally predicted only on a three- 
dimensional model. Mapping of the earth in wide spectral bands will continue to 
be a most useful meteorological tool; but, to understand the complex behavior of 
the atmosphere, a third dimension is required. 


SCIENTIFIC OBJECTIVES OF THE IRIS EXPERIMENT 
Basic Laws of Planetary Atmospheres 


Planetary atmospheres tend to approach a state of thermal and mechanical 
equilibrium, so that atmospheric phenomena may be described by these equilib- 
rium conditions and small deviations, if required. Terrestrial experience shows 
that even in strong storms the surface pressure fluctuates afew percent while tem- 
peratures between summer and winter vary less than ±10 percent. In general, 
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the atmosphere can be described on a gross scale by average or mean conditions 
combined with certain perturbations. 

The initial reaction of a hypothetical gaseous envelope surrounding a planet 
would be its adjustment to mechanical or hydrostatic equilibrium; in other words, 
pressure forces will balance gravitational forces. 

If the gas may be assumed to obey the perfect gas law, the pressure and 
density are then determined everywhere in the atmosphere. The temperatures 
are governed by the condition of thermodynamic equilibrium which requires the 
divergence of the energy flux to balance the energy sources. This net energy 
flux may be accommodated by various means of energy transport. Radiation, 
convection, heat conduction, diffusion, and even acoustical propagation of wind 
noise may participate in carrying the flux. Because many atoms and molecules 
in planetary atmospheres interact strongly with the radiation field, the radiative 
flux is very important. The radiative flux, which is always present, determines 
the thermal structure of large sections of the earth's atmosphere. When it is 
the only important mode of energy transport, a state of radiative equilibrium is 
achieved. 

Sometimes, super-adiabatic temperature gradients are required to satisfy 
the flux divergence condition. The Schwarzschild stability criterion is then 
violated, and as a consequence, convective motions occur. The slow upward and 
downward transport of gas carries the convective component of the flux. The 
zone where convection supports the radiative transport, and determines the tem- 
perature distribution, is the zone of convective equilibrium. In general, the 
earth's stratosphere may be considered in radiative equilibrium and the tropo- 
sphere in convective equilibrium. 

However, on a rotating planet the equilibrium conditions may never truly be 
reached. The difference in the net flux between polar and equatorial regions re- 
quires lateral transport of energy which leads to horizontal motions, subject to 
Coriolis forces. On earth, solar radiation is partly used to evaporate water 
vapor while precipitation and evaporation create temporary sources and sinks 
of energy in the atmosphere. The complexity is further increased by the forma- 
tion of clouds which influences the condition of thermal equilibrium strongly. 

To deal with deviations from equilibrium, the simple treatment using the 
ideal gas law, the hydrostatic and thermodynamic equilibrium conditions is in- 
sufficient. The gas law must be augmented to account for the phase changes of 
water and the formation of clouds. The Navier-Stokes equation for a compres- 
sible, viscous medium must replace the hydrostatic condition. The Navier- 
Stokes equation is a nonlinear differential equation of second order which can 
only be solved by numerical integration. 
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The integration of the Navier-Stokes equation will gain further importance 
in meteorology. At present, numerical weather-prediction techniques are based 
on the integration of greatly simplified forms of this equation. To be truly use- 
ful for weather prediction schemes, satellite data must be applicable to the 
methods of numerical forecasting. This is discussed further in the chapter on 
design considerations for this experiment. 

In the troposphere of the earth, the equilibrium conditions are of little interest. 
They are well known and only deviations from the average indicate the weather. In 
planetary atmospheres, however, the average conditions are not well known at all; 
therefore, equilibrium calculations and measurements serve very useful purposes. 


The Thermal Emission Spectrum and Planetary Exploration 


In the domain of the thermal emission spectrum of the earth and planets, a 
number of tasks can be performed. Some of these tasks such as the identifica- 
tion of gaseous constituents are very important for a study of the atmospheres of 
Venus, Mars, or Jupiter but are of little value to meteorology on earth. However, 
the known conditions on earth may serve as a proving ground for planetary ex- 
periments. Techniques as well as instrumentation may be tested. 

Information about the distribution and total amount of some constituents 
(H 2 0, e.g.) is equally important for the earth and the planets. The derivation 
of the temperature profile is again equally essential for the understanding of 
all atmospheres. The requirements on the accuracy, however, are more severe 
on earth, simply because more is known about our atmosphere than about the 
atmospheres of the other planets. The tolerances quoted in Table I give the 
order-of-magnitude of accuracy necessary to improve existing knowledge. 

On earth the surface pressure must be measured very accurately if it is to 
be of any meteorological use, because only the small deviations from average 
are of interest. For Mars, on the other hand, one would be very happy to know 
the surface pressure within 10 percent. Clearly the earth and each of the planets 
have somewhat different requirements but it will be shown later that the instru- 
mental requirements are more uniform than might be expected. 


Temperature Profile 


One objective of this experiment is to measure atmospheric as well as sur- 
face temperatures. Early investigators realized that the thermal emission spec- 
tra of stellar and planetary atmospheres contain information on the temperature 
profile. In 1956 King suggested the measurement of the limb function to determine 
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Table 1 

Parameters Influencing the Thermal Emission Spectrum 
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atmospheric temperatures on earth from a satellite and in 1959 Kaplan suggested 
a spectral scan in the 15 p. C0 2 band for the same purpose. King's method suf- 
fers from the lack of horizontal stratification in the earth atmosphere. Kaplan's 
method has been given more attention - Wark (1961), Yamamoto (1961), King 
(1963), and others - and a spectrometer was designed specifically for this pur- 
pose (Dreyfus and Hillary, 1962). 

Qualitatively the determination of the temperature profile is readily under- 
stood. The specific intensity emitted by the atmosphere originates from an 
atmospheric layer where the optical depth is about unity. Intensity is related to 
the true temperature by the Planck function (under the condition of local thermo- 
dynamic equilibrium); the atmospheric temperature can therefore be inferred 
from the intensity measurements. The resulting temperature is a weighted mean 
value for a relatively thick layer ( ~ 10 km). Using the intensity emitted in a 
spectral range of different absorption characteristics, the mean temperature of 
another layer may be obtained. In principle, then, the temperature distribution 
over a substantial portion of the atmosphere may be derived from measurements 
using a judicious selection of spectral intervals. 

Quantitatively the above argument may be understood from the solution of 
the equation of radiative transfer. 


I 


B.(T s )e- 



B (r)e Tj/ / i dr 

V' ' VfX 


( 1 ) 


For a constituent of known and uniform mixing ratio, the optical depth is a known 
function of pressure in each spectral interval. In general, r is a function of 
wavenumber, zenith angle, and the pressure and temperature corrected optical 
path. In the region of the atmospheric window, where the total optical depth to 
the surface is very small, the intensity is primarily given by the surface tem- 
perature, assuming an emissivity near unity. For opaque intervals, T is large 
and the first term in the above equation contributes nothing. If the black surface 
is considered an isothermal extension of the atmosphere, equation (1) becomes 


I 

Vfl 



B (r)e“ T ^ dr^ 


(2) 


Over narrow wavenumber intervals Av, B may be considered independent of v. 

Equation (2) is a Fredholm integral equation of the first kind and e” T is its 
kernel. The problem of finding the temperature profile requires the inversion 
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of equation (2). From the known kernel and the measured intensity, B (r) is to 
be found. It is evident from the form of equation (2) that a scan in frequency 
(Kaplan) and a scan in viewing angle (King) are formally identical problems. 

Small uncertainties in the measured intensities may result in oscillatory 
solutions of large amplitude; therefore, other criteria must be evoked to obtain 
a physically meaningful result. The problem is to apply sufficient smoothing to 
avoid instability, but not an excessive amount that would destroy the details in- 
herent in the measurements. Several methods have been discussed in the liter- 
ature on how numerical inversion can be applied to this problem — Kaplan (1961), 
Yamamoto (1961), Wark (1961), and King (1964). Flemming, Twomey and Wark 
(1964) and Dr ayson (1963) have performed error analyses and have shown that 
stringent requirements are placed on the absolute and relative accuracy of the 
intensity measurement. 

To obtain four or eight temperature levels in the atmosphere, it is sufficient 
in principle to measure the specific intensity in four or eight narrow intervals, 
strategically located so that the individual absorption coefficients provide an 
optical depth of about unity at desired pressure levels. The interferometer 
yields a continuous spectrum and, therefore, contains inherently a higher degree 
of altitude resolution than measurements in individual spectral intervals. Some 
of this information will be redundant because the absorption coefficients in many 
spectral intervals will be approximately the same. A certain degree of redun- 
dancy, on the other hand, may be used to improve the accuracy. 

The theory of the temperature sounding from satellites has been derived 
for absorbing atmospheres only. Scattering processes are usually avoided in 
calculations for reasons of simplicity. Atmospheric particles such as aerosols 
or cirrus clouds introduce difficulties which will have to be studied more 
thoroughly. 


Water Vapor and Ozone Distribution 

Suppose the analysis of the specific intensities in the 15/x band of C0 2 and 
in the atmospheric window has been performed, and the atmospheric tempera- 
tures are known down to the surface. Equation (2) may be integrated by parts 
and by a change in variables will produce: 
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= B(r 


V 



e 


t,,(T ) ^B( T , y) 

c)T 


dT 


(3) 
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Aside from B (t = 0) , equation (3) is formally identical to equation (2). The known 
kernel is now Bb/BT and one may solve for r(T). 

Numerical methods similar to those in the temperature problem can now be 
applied to a spectral region where water vapor or ozone are the major atmos- 
pheric absorbers, and their vertical distribution may then be obtained. Although 
error analyses have been performed on the temperature problem (Wark 1961, 

Dr ay son 1963), an error analysis on the determination of water-vapor distribu- 
tion has not yet been carried out. However, preliminary calculations confirm 
the feasibility of this task. King (1963) has shown a method which derives mix- 
ing ratios without the explicit determination of the temperature profile. 

Because of limited accuracy of stratospheric temperature estimates and 
greater uncertainties in the absorption characteristics, the ozone abundance and 
distribution are more difficult to determine. Earthbound and balloon-borne 
measurements of the ozone band and simultaneous stratospheric temperature 
soundings must supplement satellite -borne ozone soundings in the beginning. 

The distribution of water vapor may be derived by studying the 6.3 m band 
or the rotational bands which show sufficient strength beyond 20 m. The ozone 
distribution may best be inferred from measurements near 9.6 m* The 14 m 
ozone band is strongly blended by the 15 m C0 2 band and in addition is much 
weaker than the 9.6 ^ band. 

Absorption coefficients for water vapor and ozone are available (e.g., Goody 
1964). A verification of satellite data with simultaneous balloon soundings will 
occasionally be required to establish confidence in the conclusions. 

As in the case of temperature soundings, clouds and aerosols introduce 
difficulties in the interpretations. The seriousness of these effects must be 
studied by experimentation from high-altitude balloons. 


Surface and Cloudtop Temperatures 


After the atmospheric constituents are identified, the effect of atmospheric 
absorption can be removed from the spectrum in those spectral ranges where 
the atmosphere is not opaque. The specific intensity emitted by the surface is 
then obtained: 


I v (surface) = e v B v (T 

surface ^ 


(4) 
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In the far infrared, most natural surfaces known on earth (rocks, gravel, 
plants, snow, etc.) have a high emissivity over large spectral ranges. However, 
the crystalline structure of complex molecules create stretching and bending 
resonances. For example, quartz exhibits strong stretching bands between 800 
and 1100 cm - 1 caused by the Si-O bond; other stretching groups (Si-Si) appear 
between 600 and 800 cm -1 . The same mineral has a bending mode (Si-O-Si) at 
430 to 460 cm -1 . These resonances give a pronounced rise in the infrared re- 
flectivity which can best be observed on polished samples. The increase in re- 
flectivity is necessarily associated with a reduction in emissivity over the same 
spectral range. 

The sharp rise in reflectivity has been well known since the early days of 
infrared studies. Rubens used the selective reflectivity of quartz plates to iso- 
late narrow spectral bands in the far infrared. The reststrahlen phenomenon 
depends somewhat on the temperature of the sample and, to a much higher degree, 
on the surface roughness. Lyon (1963) has studied the emission and reflection 
spectra of many minerals in great detail. Deviations from an emissivity of unity 
in the region of the atmospheric window can cause errors in radiometric surface- 
temperature measurements. A high spectral resolution allows the identification 
of reststrahlen effects and therefore provides the means for improving the ac- 
curacy of surface-temperature measurements. This is also important for other 
radiometric, possibly scanning type, experiments especially over desert areas 
where Si0 2 reststrahlen phenomena may cause substantial errors. 

If clouds are present, the cloudtop temperature is obtained instead. By com- 
parison with the temperature measurement, the cloud altitude may be derived. 


Other Investigations 


Temperature sounding, the determination of the water vapor and ozone dis- 
tribution, and the accurate measurement of the temperature of surface materials 
do not exhaust the potentials of the IRIS experiment. A few other possibilities 
shall be mentioned. 

An accurate knowledge of the surface radiance and the atmospheric tempera- 
ture allows, in theory, a derivation of the surface pressure. As King (1963) has 
shown, a region of weak CO 2 absorption may be used for this purpose. A pre- 
liminary investigation has shown that the required accuracies cannot be met at 
the present time; however, the task is important enough to be pursued further. 

Minor atmospheric constituents such as N 2 0 (7.8; 9.56 and 17 /x) and CH 4 
(7.66/2 ) may be identified in the thermal emission spectrum. The weak absorption 
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of N 2 0 near 17 p. may be used to determine cloud altitudes in conjunction with 
the relatively transparent regions near 18 i±. 


DESIGN CONSIDERATIONS FOR THE IRIS EXPERIMENT 
Required Spectral Resolution 


The identification of gaseous constituents requires high spectral resolution. 
It would be desirable to identify individual lines in the generally very complex 
spectra. In the P and R branches of the 15 /j, C0 2 band, for example, individual 
lines have a separation of about 1.6 cm -1 and a line width of ^10 -1 cm -1 , de- 
pending on the pressure. To identify individual lines, an instrument would have 
to resolve to better than 0.1 of a wavenumber. For the earth's atmosphere, 
the identification of gases is not too important; however, the determination of 
the temperature profile requires the same order of resolving power. Unfortu- 
nately, small spacecraft instruments cannot now achieve this degree of resolu- 
tion in the infrared. 

The next best approach is to resolve the band contour as well as possible. 
Then a spectral resolution which is small compared to the total bandwidth is 
sufficient to identify a particular constituent. Similarly, for the measurement 
of the temperature profile, groups of lines may be selected. Kaplan (1959), Wark 
(1961) and others have shown that a five-wavenumber interval is sufficient. 
Smaller intervals do not yield further improvement unless one reaches 0.1 of a 
wavenumber or less and resolves individual lines. A five -wavenumber interval 
is within the feasible range for satellite instrumentation. 

In summary, the goal of five-wavenumber resolution adopted for this exper- 
iment is a compromise between the desirable capability of 0.1 wavenumber and 
instrumental limitations. However, five-wavenumber resolution is a substantial 
improvement over space-borne techniques presently available and is sufficient 
to answer many pressing questions. 


Required Spectral Range 


The desired spectral range will be derived for the atmospheres of the earth 
as well as for other planets. However, only the region of the thermal emission 
spectrum where the reflected solar radiation is negligible will be considered. 
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• Earth 


As pointed out earlier, a very interesting task will be the measurement 
of the temperature profile. Consequently the absorption spectrum of a uniformly 
mixed gas such as C0 2 should be included (15 p). After deducing the tempera- 
ture profile, the abundance and distribution of gases of variable mixing such as 
H 2 0 and O3 may be obtained. Therefore, water-vapor and ozone absorption 
regions are of great interest (6.3 p, 9.6 and the rotational bands of water 
vapor). The surface or cloudtop temperatures may be measured in the atmos- 
pheric window between 10.5 and 11.5 /x, and this region should therefore be in- 
cluded in the desired range. Furthermore, if the surface and atmospheric tem- 
peratures and the surface emissivity are known, a region of intermediate C0 2 
absorption may ultimately allow a direct determination of the surface pressure. 
However, present measurements are far from the required precision. 

Except for small gaps, the whole thermal-emission spectrum yields signifi- 
cant information. Since the proposed interferometric method allows the cover- 
age of a wide spectral region, and since the exclusion of certain gaps is more 
difficult than transmitting the whole spectrum, a record of the emission spectrum 
from 2000 cm -1 (5/x) to 500 cm -1 (20/u) is proposed. It would be desirable to 
cover the strong CO 2 bands near 4.2 /x , also, but the signal-to-noise ratio of the 
present instrument becomes very poor in this region. A response beyond 20 /x 
would also be desirable; the rotational water-vapor band would be superior over 
the 6.3 fx region from a signal-to-noise point of view, but other instrumental con- 
siderations (materials for beam splitter) impose limits there. The spectral re- 
gion of the instrument will hopefully be extended to 30/x in the future. 


• Venus 

The abundance of C0 2 on Venus makes the measurement of the C0 2 bands 
very important including the weaker bands of C0 2 near 9.4 and 10.4 p. Water- 
vapor absorption in the 6.3 /x band would be sufficient to confirm present estimates 
on the water vapor content above the clouds. Today, nothing is known about the 
thermal emission of Venus between 12 and lOOp where 75 percent of the thermal 
energy lies. Besides the search for known or expected phenomena, a search for 
unexpected spectral features is most important. Therefore, a wide spectral range 
should be the design goal for a Venus fly-by experiment. 


• Mars 

Similar arguments hold for a possible exploration of Mars. In addition 
to the identification of atmospheric constituents, " reststrahlen" phenomena in 


10 



surface materials may allow their identification. Again, a search for unexpected 
features requires a wide spectral range (Hanel and Chaney, 1964). 


• Jupiter 

The atmosphere of Jupiter contains large quantities of hydrogen and he- 
lium, but neither molecule shows resonances at moderate pressures in the far 
infrared region. However, other constituents such as NH 3 and CH 4 have strong 
molecular absorption bands in the 5 to 30 p range. Again, a wide spectral range 
is most desirable. 

In summary, a wide spectral range of the instrument is most desirable for 
application to the earth's atmosphere for meteorological purposes and for plane- 
tary exploration from fly-by vehicles. As mentioned before, the design goal en- 
compasses a spectral range from 5 p. (2000 cm -1 ) to 20p (500 cm -1 ) with a pos- 
sible extension to longer wavelength. 

The wide spectral range does create challenging problems in the design of 
the instrument, primarily in two areas: One concerns the beam splitter, which 
has to operate properly over the wide range. The other problem is in the large 
dynamic range of the signal which must be digested by the data transmission 
channel. As discussed later, both problems have been solved satisfactorily. 


Required Spatial Resolution in the Horizontal and Vertical 

Weather prediction by numerical integration of the Navier-Stokes equations, 
or simplified forms of it, requires the knowledge of certain initial conditions. 
Initial conditions reflect the state of the atmosphere at a particular time and 
consist ideally of the wind velocities, the temperature profile, the water vapor 
content, and the surface pressure. For short time forecasting, on the order of 
a few days, atmospheric heating processes are less important than the afore- 
mentioned parameters (Mintz 1963). With the exception of the surface temper- 
ature, and those rare occasions when the wind velocities can be inferred from 
cloud patterns or from the size of the sun's image reflected from water, neither 
of the necessary parameters is measured by meteorological satellites today. 
Wind fields and surface pressure are still beyond the realistic capabilities of 
remote measurements, but the temperature profile and the water vapor distri- 
bution are within the potential of this experiment. 

For longer range forecasting, the time -dependent heating or cooling of the 
atmosphere also becomes important. Surface temperatures and cloudiness are 
main factors. Both these parameters, in addition to the total radiative heat loss, 
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are already under satellite observation or can be inferred from the results of 
such observation. 

Present techniques for numerical forecasting are still relatively crude. 

The finite storage capability of today's computers limits the number of atmos- 
pheric volume elements which can be considered simultaneously. As discussed 
by Mintz (1963), the lowest useful resolution uses a two-layer model of the 
atmosphere with a horizontal dimension of 1000 km per element. This type of 
element is still small compared to the characteristic wavelength of cyclonic 
motions. Surprisingly, numerical integration of the so-called primitive equa- 
tions yields very good results using a two- or four-layer atmospheric model 
with a horizontal element size of only 800 km (Mintz, 1964). 

These considerations are useful for comparison with the minimum spatial 
resolution requirement for the interferometer experiment. Spectral resolution 
and wavelength range limit the maximum field-of-view of the interferometer to 
about 8 degrees. This is discussed in detail later. At the nominal altitude of 
Nimbus (1100 km), this field-of-view corresponds to a circular area of about 
160 km diameter. The horizontal resolution of this instrument is therefore much 
better than the element sizes required for numerical forecasting. The interfer- 
ometer will be used without additional optical elements on Nimbus. Other ap- 
plications such as planetary fly-by missions from greater distances, or observa- 
tions from synchronous altitudes on earth, require a telescope to take full 
advantage of the light gathering power of the interferometer and simultaneously 
restrict the overall field of view to less than 8 degrees. 

The number of vertical layers for which the temperature and water vapor 
content should be derived must again be determined by the number of layers that 
the numerical models of atmospheric behavior can digest. From the present 
state-of-the-art of numerical forecasting, it is concluded that four or perhaps 
eight layers distributed equally over pressure, if possible, should be ample. 


Dynamic Range 


The IRIS instrument, like all radiometric devices using thermal detectors, 
responds to the difference between the radiance of the target and its own radiant 
emission. Three principal modes of operations are possible: The instrument 
temperature can be chosen either entirely above or below the target temperatures, 
or near the center of the dynamic range. Operation at either temperature extreme 
avoids ambiguities in the interpretation of the spectra at the expense of the re- 
quired dynamic range. Since outer space will be used as a convenient almost zero- 
degree calibration target, operation below target temperatures is impossible. 
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Operation above the maximum of expected blackbody temperatures (320°K in 
window over desert) is undesirable, because of the wide dynamic range required 
and because of noise considerations in the detector. Therefore, an operating 
temperature at 250 °K was selected. Under these conditions the emitted radia- 
tion during calibration approximately equals the received net radiation from the 
warmest target (300°K). It is sufficient to design for the maximum average 
blackbody temperature which may be encountered, rather than for the peak, since 
phase reversals are permitted. However, the insensitivity of the power spectrum 
to the relative "phase” of the individual frequency components in the interfero- 
gram leads to an ambiguity. The power spectrum was selected rather than a pure 
cosine transformation for this very reason. We do not wish to depend on the phase 
characteristic of the system and the lack of dispersion in the beam splitter. On 
the contrary, a moderate dispersion is desired since it reduces the high central 
peak in the interferogram. 

In general, the ambiguity can be resolved by an a priori knowledge of the 
shape of the spectrum; however, in automatic modes of data processing, a more 
definitive criterion is required. This criterion exists in the difference in "phase" 
of two adjacent frequency components; if they appear "in phase" near the center 
region of the interferogram, the corresponding blackbody temperatures are both 
warmer or both cooler than the instrument temperature. An accurate knowledge 
of the position of the central peak is not required. Because certain parts of the 
spectrum such as parts of the 15 / j . C0 2 band will always be cooler than 250°K, 
the ambiguity is resolved. The relative "phase" of the individual frequency com- 
ponents is obtained from the sign of the sine and cosine amplitudes which are 
squared and summed to yield one data point in the power spectrum. 


THE INTERFEROMETER SPECTROMETER 

The instrument, basically a Twyman-Green modification of the Michelson 
interferometer (Figure 1), departs from the classical laboratory designs in 
several areas: 

• The drive mirror is suspended by a parallel leaf-spring system. 

• The mirror assembly is driven by a moving coil in a strong magnetic 
field; a second coil generates a voltage proportional to the mirror speed 
which is used to control the mirror speed by feedback. 

• The mirror motion is monitored by an auxiliary interferometer to as- 
sure proper sampling times. 

• The compensating plate is moved to a position adjacent to the beam splitter. 
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Figure 1 - Interferometer Cube, Mirror Drive Assembly, 
and Detector Assembly: Schematic Diagram 


This particular configuration of the interferometer was chosen because of 
its inherent compactness, large light-gathering power, and design simplicity. 

Radiation arriving within the field-of-view is divided by the beamsplitter 
into two components which recombine and interfere after reflections on the 
fixed and moving mirrors respectively. The recombined beam, focused on the 
detector, is modulated by the motion of the scan mirror proportional to the speed 
of the mirror and the wavenumber of the incident radiation. The output signal 
from the detector, called the interferogram, is the Fourier transformation of 
the spectrum of the incident radiation. The signal is amplified, digitized, and 
transmitted over the spacecraft telemetering system. On the ground, the spec- 
trum is reconstructed by applying the inverse transformation. The original 
spectrum can be precisely reconstructed only within certain limits of accuracy 
of spectral range and of spectral resolution. 
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Theory of Fourier Spectroscopy 


Interference spectroscopy, first suggested by Michelson after the develop- 
ment of his original interferometer, was not pursued because of extreme com- 
putational difficulties, and further development was forced to wait upon the 
advent of electronics and high-speed digital computers. In the last decade, in- 
terference spectroscopy has made great advances; work has been reported by 
Fellgett (1958), Gebbie (1959), Mertz (1956), Connes (1961), Strong and Vanesse 
(1959), Loewenstein (1963), and many others. Excellent summaries may be 
found in articles by Jacquinot (1960) and by Connes (1961). 

The source of radiation for the interferometer will be an extended target 
which completely fills the field-of-view with radiation passing through the in- 
strument without collimation. The output fringe pattern of the interferometer 
is focused on the detector, which is masked to eliminate all but the central fringe. 
The output of the detector is amplified and recorded as a function of path differ- 
ence, d, between the fixed and movable mirror. 

The relationship between the effective path difference, the order of inter- 
ference, and the radiation wavelength is 

m\. = 2d cos 6 


where 

d is the path difference between the two interferometer paths (cm) 
m is the order of interference 
A. is the wavelength (cm) 

6 is the angle of an oblique ray through the interferometer. 


This expression shows that, for a given monochromatic signal of wavelength, 
A., the circular fringe corresponding to a given order m expands away from the 
center of the pattern as the separation of the mirrors is increased. The detector 
is mounted behind the mask so that only the central fringe is detected; each spec- 
tral element is thereby modulated sinusoidally at a rate inversely proportional to 
its wavelength. As the source has an appreciable range of spectral elements, the 
signal detected will be the superposition of many modulated intensities multiplied 
by a parameter which is the "responsivity" of the system; i.e., 




(X) 




d V 


K N 


COS 2TTUxdv 


15 



where N v is the spectral radiance of the source (assumed to be the same at all 
parts of the field-of-view) 


v is the wavenumber, cm 1 (v - 1/k) 


x is the instantaneous value of the effective path difference of the two 
beams of the interferometer (cm) 

K„ is a parameter which expresses the effective spectral responsivity 
of the interferometer, including the characteristics of the optics, 
the detector, and the electronic amplifiers 


The interferogram is then taken to be 


I 


(*) 




K v N v cos 2rrvxdv , 


which is the Fourier cosine transform of K v * N v . The product K„ • can be 
obtained uniquely from the measured interferogram, taking the inverse trans- 
form by digital computation. The relationship between N v and the product 
K„ ■ N„ will be established by laboratory calibrations; the validity of the labo- 
ratory calibrations is established by in-flight calibrations. 

The true spectral distribution • N y is found only when the interfero- 
gram is recorded for an infinite range of x . For a finite range of x, each spec- 
tral component is spread out into a shape known as the instrumental profile. 
This profile also depends on the obliquity of rays through the interferometer, 
and thus on the solid angle of the field-of-view. The instrument profile may be 
modified as desired during computation of the transform; this process, called 
apodization, usually is done to remove secondary maxima at the expense of 
resolution. 

Including all these effects, the expression for the interferogram becomes 

r 

I = N v ■ s(v,x) • cos [ 2rrvx + <fc( v, x) ] , 

Jo 
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where 


S(v,x) - s (y ,x) exp [ i0(v ,x) ] 

is the Fourier transform of the instrument profile (including the parameter K„ 
mentioned previously). 

The spectral distribution is again obtained by computing the Fourier trans- 
form of I (x) . Reconstruction of the spectral distribution of the source N v is 
obviously an approximation of the actual spectral distribution. This approximate 
solution (i.e., the transform of the interferogram) is the convolution of the actual 
spectral distribution with the instrument profile; therefore, it is mandatory to 
determine the form of the instrument profile. This can be done experimentally, 
because the instrument profile is the apparent spectral response of the instru- 
ment to a monochromatic source. This determination of the instrument profile 
is an important part of the laboratory calibrations. 

Two factors which influence the instrument profile and reduce the quality 
of the data are the limitation arising from the finite path of the movable mirror, 
and the effect of oblique rays through the interferometer. The resolution limit 
imposed by the finite-path limitation is 


R = vL or 



(assuming linear apodization to remove secondary maxima in the instrument 
function). The resolution limit imposed by oblique-ray effects is 


RQ = 277, 


v 0. 

or Av = 1 “ ’ 


In the above, R is the resolution v/Av, L is the maximum path difference on 
either side of zero path difference, and 0 . is the solid-angle field -of-view of the 
interferometer. 

• Advantages of the Interferometer 

An analysis by Jacquinot indicates that, where a source whose spectrum 
is defined by a radiance N , having M spectral elements of interest, is to be 
analyzed with a resolution R in a time T with an optical system which is 
detector-noise-limited (as presently in the infrared), then the signal-to-noise 
ratio is increased by maximizing the time t of observation of each spectral 
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element, the admissible solid angle Q, and the transmission factor e v (the opti- 
cal efficiency). 

An interferometer fulfills these requirements to a high degree. The time 
of observation is t = T(m/M) ,, where m is the number of spectral elements 
studied simultaneously. For an interferometer, m equals M , and m/M has the 
maximum value 1. The luminous energy received by the detector is N v • A • fi . 
The quantities A and Q can both be made large for an interferometer, as long 
as the resolution limitation RQ = 2n is not exceeded (for comparison, the same 
resolution limit for a grating spectrometer is Rfi < 0.1). The transmission 
factor « is the only area where the superiority of the interferometer may be 
questioned for a given application. The reflectivity and transmissivity of the 
beam splitter is of primary importance, and it may be difficult to obtain high 
efficiency from mechanically suitable materials in certain portions of the infra- 
red spectrum. 


Mechanical and Electromechanical Design of Breadboard Instrument 


The instrument for the IRIS experiment consists of: 

• interferometer (cube, detector, mirror-drive system, preamplifiers) 

• 45-degree mirror rotating to reflect outer space and a blackbody during 
in-flight calibration and to provide image-motion compensation 

• electronics compartment, containing all electronic circuits except for 
the preamplifiers in the cube and electronics located in calibration mir- 
ror assembly 

• mechanical structure to receive the interferometer, the calibration 
target, the calibration mirror, and to provide mechanical strength and 
thermal protection 

The electronics compartment, located separately, is connected to the interfer- 
ometer assembly by a cable. 

The mechanical design of the interferometer is modular. The drive unit, 
the fixed mirror assembly, the detector assembly, and the monochromatic ref- 
erence source-detector combination are integral components attached to the 
optical cube, which houses the beam splitter and provides mechanical strength 
and thermal conduction to all components. 
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The scan mirror is one of the most critical items in the design because it 
must be moved through a distance of 2 mm with a total angular rotation of less 
than 2 seconds of arc. 

The use of parallel leaf springs to support the mirror was chosen in prefer- 
ence to slides on ways, or pivot-point mounts, since it avoids frictional surfaces 
undesirable in space applications. The driving force is generated by an electrical 
current in a moving coil in a strong magnetic field. As Figure 1 shows, the mag- 
netic gap is inside the dual-magnet design, thus minimizing residual magnetic 
fields on the outside. 

Lack of uniformity in drive-mirror speed normally produces catastrophic 
effects on the wavenumber calibration of interferometers. Two methods used 
to eliminate virtually all errors and uncertainties in the mirror motion, and in 
the subsequent wavelength calibration, are based on the use of an auxiliary inter- 
ferometer and a feedback control of the mirror drive: 

• A monochromatic source (0.5852 micron) feeds the interferometer 
and generates a sinewave at an auxiliary detector. The frequency of 
the sinewave is directly proportional to the speed of the mirror. 

Sampling times of the interferogram are derived from this reference 
frequency. 

• The mirror velocity is maintained constant to ±1 percent by energizing 
an electromechanical drive coil from an error amplifier that compares 
a reference voltage to a voltage proportional to the mirror velocity. 

The latter voltage is derived from a pickup coil attached to the mirror 
mount that moves through a constant magnetic field. This method also 
provides electromechanical damping in the neutral position of the mir- 
ror (with zero reference voltage) which is very important during the 
launch phase. 


Optical Design of the Breadboard Instrument 


The optical design of the interferometer is conventional. The flat mirrors 
and the detector determine the limiting apertures. The condition of optimum 
solid angle RO = 2n is fulfilled at the highest wavenumber to be analyzed 
(2000 cm -1 , 5 >u. ). The optimum angle is then 8 degrees (±4 degrees) for a re- 
solving power of 400 (5 cm -1 at 2000 cm -1 ). 

The fixed mirror and the scan mirror are made of gold-coated quartz pol- 
ished to a flatness better than 0.1^. Gold was selected because of its stability 
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and its high reflectivity in the infrared. The two-axis mirror mount is machined 
from a single piece of metal. A differential screw which produces the fine ad- 
justment required also anchors each adjustable plane at three points. 

The condensing optic consists of a KBr lens, combined with a gold-coated 
reflective cone which increases the effective detector area. Dimensions are 
chosen to provide the 8 -degree acceptable angle. The f-number of this system 
will approach the ideal limit of 0.5. 

The bolometer in this design is essentially an integral part of the condensing 
optic. The detector is a 3.5-mm square masked to a circle of 3.5-mm diameter, 
having a time constant of 1 millisecond. 

The beam splitter is the most critical item in the optical train. All the rays 
make two passes through the beam splitter; hence, all the losses arising in this 
unit are doubled. Although metallic beam splitters are preferable from the stand- 
point of uniformity with wavelength, they have the serious disadvantage that ab- 
sorption losses are high. Dielectric beam splitters cannot be made uniform with 
wavelength, but over the wavelength span of interest the efficiency will be at least 
equal to that of a metallic splitter. The dielectric material selected is germa- 
nium, which has the additional advantage of being highly reflective for wave- 
lengths shorter than 2/j., so that most of the reflected visible radiation will not 
be collected. 

The substrate material used in the breadboard instrument is KBr. Other 
materials currently being investigated for their suitability as beam splitters 
include IRTRAN-4, KRS-5, Csl, and plastic pellicles. 

The radiation source for the monochromatic reference signal is the 0.5852 p- 
neon line, which is quite strong; in conjunction with a suitable interference filter, 
the signal is suitably monochromatic. Neon was selected because lamps are 
available in small sizes and very little power is required. The focusing and 
condensing optics are all glass of conventional and straightforward design. A 
silicon photovoltaic cell serves as a detector. 


Electronic Design of the Breadboard Instrument 


The electronics for the IRIS instrument are designed to perform the func- 
tions necessary to control the interferometer and process the data presentation 
to the spacecraft. The design consists primarily of digital electronics, because 
of the inherent reliability of these systems. The digital circuits are constructed 
of Texas Instrument modules. This subsystem is a self-contained instrument 
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requiring a minimum of external control for the operation of the experiment. A 
reliability survey will be conducted as one of the system evaluation parameters. 

The electronic circuitry for the interferometer is designed to operate as 
independently of the spacecraft as possible. The spacecraft supplies the power 
and command signals necessary to operate the experiment; the electronics sub- 
system of the experiment processes the interferometer data and presents them 
to the spacecraft in a digital format. The spacecraft supplies a single command 
at the beginning of each interferometer frame, but the necessary sequence of 
events to complete the data frame is generated by the interferometer system. 

Functions of the interferometer electronics (Figure 2) are: 

• To amplify the thermistor-bolometer output signal 

• To sample the bolometer signal and present it to the spacecraft telem- 
etry system in serial digital form 

• To provide the necessary control to drive the moving mirror at constant 
speed (±1 percent) 

• To provide a calibration sequence and controls for image-motion 
compensation 

• Convert the -24. 5v spacecraft power into voltages required to operate 
the interferometer 

The data channel receives as an input the 0 to 1.0-volt peak analog signal 
from the thermistor bolometer preamplifier. This electrical signal is band- 
limited between 10 cps and 100 cps in an electrical filter network, 12 db/decade. 

The signal amplifier has a voltage gain of 36 db. The band -limited electrical 
signal is coupled through the switched attenuator (0, - 20 db) into the signal am- 
plifier. The state of the switched attenuator is controlled by a comparator cir- 
cuit that senses the level at the output of the signal amplifier, as shown in Figure 2. 
The attenuator has two positions, controlled by the level of the preamplifier output 
voltage: 0-db attenuation and 20-db attenuation. The preamplifier output fed into 
a comparator generates an output to switch in the 20-db attenuator when the pre- 
amplifier output exceeds ±3.2 volts, and to switch out the attenuator when the 
signal falls below ±0.32 volts. The gain change is necessary to accommodate 
the extremely large dynamic range (1500:1) of the input signal from the bolom- 
eter preamplifier. The output of the bolometer preamplifier is maintained be- 
tween 0 and -6.4 volts, a value compatible with the analog-to -digital converter. 
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ANALOG DATA CHANNEL 
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The function of the data control channel is to: 


• Provide the necessary encode pulses to the analog-to -digital converter 
at proper mirror positions 

• Ensure that data points are not transmitted until the mirror velocity has 
stabilized 

• Allow precisely 3408 data points for each data frame 

• Allow 64 additional pulses to cycle through the 16 housekeeping telemetry 
points four times 

The encode pulses to the analog-to-digital converter are derived from the 
auxiliary interference of the neon reference. This ensures that the data samples 
are taken at specific mirror positions rather than at equal intervals of time. 

This precaution is taken to remove the effect of small mirror-velocity variations. 

The count-to-128 logical output inhibits data flow from the experiment until 
the mirror velocity has stabilized. After this time, the 312-1/2 cps neon inter- 
ference signals are coupled to the analog-to-digital converter as encode pulses. 

The logical output from the count -to-3408 switches the operation of the in- 
strument from the infrared data to the housekeeping data. The 16 housekeeping 
channels are cycled through four times, and the instrument is switched off. 

The telemetry output unit receives the digital data from the analog converter 
and electronics in the form of a 12 -bit word. The word consists of the 8 bits from 
the analog-to-digital converter and a parity bit, a single bit indicating the state of 
the 20-db attenuator, and two bits used for word synchronization. These data are 
transferred to the spacecraft tape recorder by a clock that is phase-locked to the 
neon interference signal. 

The function of the mirror-control circuitry is to drive the Michelson mirror 
at a constant velocity for the duration of one experiment. The design goal for the 
mirror-velocity control servomechanism is for a velocity of 0.185 mm/sec ±1 
percent for a displacement of about 2 mm. The velocity is kept constant by con- 
trolling the current in an electromechanical drive coil attached to the moving 
mirror. This current is derived from an operational amplifier coupled to a 
reference voltage and a voltage proportional to mirror velocity. This voltage 
proportional to velocity is generated in a pickup coil attached to the mirror 
mount that moves through a constant magnetic field. 
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The 45-degree scan mirror compensates for image motion to ensure that 
the same scene is viewed during the 11 -sec interval of each interferogram. 

In addition to providing motion compensation, this scan mechanism provides 
calibration signals. During calibration the scan mirror is rotated to bring a 
thermally controlled (300°K) blackbody into the field-of-view of the instrument. 
After obtaining the calibration interferogram, the scan mirror rotates 90-degrees 
to bring outer space (0°K) into the field-of-view and a second calibration interfer- 
ogram is recorded. This constitutes calibration for both extremes of the instru- 
ment. The calibration sequence, which will be initiated after every 14 interfer- 
ograms (every 16th frame-initiate command), is inhibited by a signal from the 
solar panel position indicator to guard against damaging the sensor by projecting 
the solar image on the detector. 

The power supply is a standard type of dc-to-dc converter. The raw voltage, 
necessary to operate the instrument are regulated to ±5 percent, sufficient for 
proper operation. In addition two reference voltages are regulated to ±0.1 percent. 
Outputs from the power supply are: 

• ±300v ± 0.1 percent thermistor bolometer bias supply 

• 12v ± 0.1 percent reference voltage - mirror drive 

• ±18, ±12, ±6, ±4v at 2.0 percent regulation 


Signal-to-Noise Ratio 

The peak signal falling on the detector (for zero path difference) is 

P = -jN • A • fi • e . 

where 

P = radiant power (watts) 

N = target radiance (watt cm” 2 ster -1 ) 

A = effective aperture (cm 2 ) 
fi = solid angle (ster) 
e = average optical efficiency 
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The factor 1/2 accounts for the 50-percent reflection losses inherent in this 
type of interferometer. The noise at the detector is given by the noise-equivalent 
power (NEP) of a thermistor bolometer, which is approximately 


NEP = 1.6 



where 


a = area of detector (cm 2 ) 

Af = noise bandwidth (cps) 
r = detector time constant (sec) 

Additional noise will be contributed by the compensating thermistor (/2), 
noise in the bias supply (/2), and preamplifier noise (/2). For operation at 
250°K the NEP is lower; however, this improvement will not be considered here. 

The peak signal-to-RMS noise ratio in the interferogram, for zero path 
displacement, is then 


S_ 

N 


P 

NEP 


10 9 NAfie 




aAf 


The signal power in the spectral interval Ay will be approximately equal to 
the peak power P times the ratio N V A^ e^/Ne . In addition, a signal -to-noise im- 
provement factor proportional to Vn/2 (n = number of samples per interferogram) 
is obtained; therefore, 



10 9 \A v e v A 


n 


yAE. 

f 2aAf 


The square root of n/2 is used instead of /n", because about half of the samples 
are redundant. 

Assuming reasonable numerical values for some parameters, 
effective aperture A =10 cm 2 

solid angle a = 1.57 • 10’ 2 (8°) 
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optical efficiency 
detector time const, r 
detector area a 
number of samples n 
spectral element Ai/ 
noise equivalent bandwidth Af ~ 
the peak signal -to-noise ratio becomes 


= 0.2 
= 1 ms 
= .1225 cm 2 
= 3408 
= 5 cm' 1 
136 cps 


— - 2. 5 • 10 s N, 
N 


and for a spectral interval 



10 7 N v . 


The instrument temperature (250°K) was chosen to give approximately equal 
signals for the maximum blackbody temperature to be measured (300°K), and for 
outer space (0°K), which is used as a check of calibration. The maximum radi- 
ance between 500 and 2000 cm’ 1 is about 6.2 • 10’ 3 watt cm" 2 ster 1 (300°K); 
therefore, the maximum signal -to-noise ratio is about 1500. 

Within the spectral intervals, the signal -to -noise ratios depend largely on 
spectral radiance. A plot of N y for various blackbody temperatures is shown in 
Figures 3 and 4; in both cases, the instrument is at 250°K. However, before the 
final signal-to-noise ratio can be discussed realistically, errors in the coding 
process must be taken into account. 

The most efficient way to encode is to set one coding step equal to the rms 
value of the noise in the interferogram. Higher coding accuracy yields insig- 
nificant improvement; lower coding accuracy is undesirable. If encoding is per- 
formed to 1 part in 1500 of the maximum possible amplitude in the interferogram, 
the uncertainties introduced by the coding process will equal errors caused by 
detector noise. The final signal-to-noise value, including the coding error, 
becomes 
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Figure 3— Difference in Spectral Radiance of a Target Colder Than the Interferometer 


— = 1.7 • 10 s N and ® = 3.5 • 10 7 N V . 

N \N/Ay 

A (S/N)a„ equal to unity is obtained for N v = 3 * 10 8 . This value is indicated 
in Figures 3 and 4. A temperature difference of 2°K should therefore be detect- 
able below 1650 cm’ 1 (above 6 microns); at 1850 cm' 1 (5.4 microns), 4°K cor- 
responds to the noise, and at 2000 (5 microns) only about 8°K. 


27 




Figure 4-Difference in Spectral Radiance of a Target Warmer Than the Interferometer 


These numbers were derived for 5 cm" 1 resolution. Clearly, numerical 
smoothing techniques can be applied and resolving power can be traded against 
accuracy. This operation is done on the ground after the interferogram has 
been received, so that considerations of this nature do not influence the design 
of the instrument. 

It remains to be shown how a signal-to-noise ratio of 1500 can be obtained in 
the coding process. The large central peak in the middle of most interfero grams 
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is very pronounced; in most cases, the amplitudes in the rest of the interfero- 
gram are considerably smaller. An 8-bit word is used to encode on a scale of 
plus and minus 128 levels, or a total of 256. One level is equal to the rms value 
of the noise. Larger amplitudes are divided by 10 electronically and then en- 
coded. The position of the divide-by-10 circuit is indicated by a ninth bit. This 
effective dynamic compression maintains full accuracy for most points in the 
interferogram, but only one -tenth of the full accuracy for all points exceeding 
a tenth of the maximum possible amplitude. However, the one-tenth value in 
amplitude is exceeded for only a very small number of data points, compared to 
the total of 3408 points. The slightly reduced accuracy of a few data points can 
be shown to have negligible influence on the reproducibility of the spectra. 


Thermal Design 


Special design considerations are required to maintain an operating temper- 
ature of 250°K. The Nimbus spacecraft temperature is maintained at about 295°K 
by thermostat. In principle, it would be sufficient to cool only the detector, but to 
minimize the effect of stray radiation, it is advantageous to maintain a uniform 
temperature in the whole interferometer cube assembly. This requires good 
thermal conduction between all elements within the cube assembly while simul- 
taneously isolating the assembly from the spacecraft and from all surfaces except 
one side which never faces the sun and the entrance aperture, of course. The 
exposed side will have a high emissivity beyond 5p and will therefore radiate 
about 2 • 10' 2 watt cm' 2 into outer space. 

Considering a total radiating area of 12 x 25 cm 2 , and allowing for small 
obscurations in the field-of-view, the total heat loss amounts to about 5 watts. 

The loss is compensated by the electrical power dissipated in the interferom- 
eter circuits located in the cube (~1 watt), by thermal energy conducted from 
the spacecraft through the mounting structure (~ 3 watts). The remaining dif- 
ference of approximately 1 watt is generated by a small electrical heater which 
thermostatically controls the instrument to 250°K ± 2°K. 


Image Motion Compensation 

The interferometer is rigidly mounted to the spacecraft. Within the 11 
seconds required for one interferogram, the orbital motion of Nimbus corre- 
sponds to about 80 km, or one-half the field-of-view. A change in radiance from 
the earth caused by the appearance or disappearance of clouds within the field- 
of-view would increase the difficulties of interpreting the spectra. The instru- 
ment registers some time-average of the spectra, but the spectral resolution is 
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also reduced. Sudden changes in the radiance at one wavelength will generate an 
amplitude modulation on the corresponding frequency in the interferogram. In 
the analysis, sidebands are generated which interfere with the signals from 
neighboring wavenumbers, thus reducing the overall accuracy as well as the 
spectral resolution. 

The image-motion compensation effectively eliminates the smear and as- 
sures that the scenery does not change during the recording of one interfero- 
gram. This is accomplished by a slow rotation of the 45-degree calibration 
mirror. At the beginning of each interferogram, the mirror images an area 
2 degrees ahead of the satellite subpoint into the interferometer. The same 
area is then kept under observation for 11 seconds. At the end of this period, 
the mirror points 2 degrees behind the satellite and advances to the starting 
position within the 5-second time interval between interferograms. 


Auxiliary Measurements 


A number of auxiliary measurements will be performed to facilitate better 
interpretation of the interferograms and to validate the operation of the instru- 
ment. 


• Temperature Sensors 

Various temperatures will be measured in the optical cube (the temper- 
ature of the bolometer mounting, the drive mirror, the beam splitter mounting, 
and a point near the entrance aperture) with an accuracy of 0.5°C. These tem- 
perature measurements will aid in the proper interpretation of the interfero- 
grams. The blackbody used for calibration purposes will be monitored by two 
bead thermistors to provide redundancy in this measurement. 


• Voltage Measurements 

A number of important voltages and execution commands will be moni- 
tored to certify proper operation of the instrument. The measured voltages 
include the bolometer bias voltage, the main power in the optical cube (-12 and 
-18 volts), and other strategically important voltages in the electronics com- 
partment. A total of 15 auxiliary measurements will be performed. 
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DATA REDUCTION 


The task of data reduction is more complicated for this experiment than for 
many others. A relatively high number of bits (almost 42,000 per interferogram) 
must be received, formatted, and presented to a digital computer for further proc- 
essing and analysis. The purpose of this analysis program is to produce a power 
spectrum which is obtained from the Fourier transformation of the interfero- 
gram. This is necessary because the interferogram itself is not subject to de- 
tailed scientific analysis, and dependence upon the phases of the spectral com- 
ponents present is undesirable. 

To test this analysis program, a synthesis program was created which simu- 
lates the performance of the instrument in detail, including the effects of obliquity 
of rays (caused by the finiteness of the optical solid angle), detector noise, and 
quantizing errors. 

Various instrumentation design parameters were determined or confirmed 
with the aid of these synthesis and analysis computer programs. A description 
of the programs and some of the results of this study have been reported by 
Byrne (1965). 

A sample, simulated interferogram is shown in Figure 5; Figures 6 to 9 
show the transformed spectral profile for each of four cases. In every case the 
same source spectral profile was used (Wark, Alishouse and Yamamoto 1963, 
Atmosphere "G", Z G = -6371); this is the spectrum for a clear warm atmosphere. 
This input spectrum is indicated on each spectral plot by the solid line segments 
shown, and the transformed spectrum is given by the "asterisked" line segments. 
Wherever noise was applied, the peak signal-to-noise (RMS) ratio was taken to 
be 1270:1, and except for Figure 9, each interferogram was computed for 3407 
equally-spaced samples over a 10-second interval. The interferogram corre- 
sponding to the transformed spectrum of Figure 9 was taken for 1707 samples 
over the 10-second interval. 

Figure 5 is the interferogram with all effects included, and Figure 6 is the 
corresponding transformed spectrum. 

Figure 7 shows the effect on the transformed spectral profile of the finite 
solid angle and of quantization only; Figure 8, the effects of finite solid angle and 
detector noise only. Figure 9 includes all effects but, as noted above, the slower 
interferogram sampling rate was used. 

In summary, results from the use of the two computer programs indicate 
that a sampling rate for I RI S of approximately four times the highest frequency 
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Response of the Interferometer to a Clear Warm Atmosphere was Calculated on a Digital Computer. Figure includes 

Effects of Finite Solid Angle, Detector, and Quantizing Noise 

Figure 5-Simulated Interferogram 
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Figure 6- Power Spectrum Computed from Figure 7 - Power Spectrum of Simulated 
Simulated Interferogram Interferogram, Including Only the Effects 

of Finite Solid Angle and of Quantization 
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Figure 8— Power Spectrum of Simulated Figure 9 - Power Spectrum of Simulated 
Interferogram, Including Only the Effects Interferogram Derived from Only 1707 
of Finite Solid Angle and of Detector Samples 




present is desirable. The effect of quantizing on the transformed spectral pro- 
file is, as anticipated, about the same as the effect of the mean expected signal - 
to-noise ratio at the detector. The gain switch, or the divide -by -ten, produces 
a negligible effect. The effect of the finite optical solid angle is less than that 
caused by digitization or noise at the detector. 


DESIGN PARAMETERS FOR INFRARED INTERFEROMETER 


Spectral range 

Spectral resolution 
Travel of drive mirror 
Diameter of effective aperture 
Detector 

Detector time constant 
Duration of interferogram 
Electrical frequencies at detector 
Sampling 

Sampling rate 

Samples per interferogram 

Housekeeping information 

Samples per frame 
Bits per word 
Bits per frame 


2000-500 cm' 1 (5 - 20^) 

(possible extension to 300 cm' 1 (33 (j.) 

5 cm" 1 over total range 

0.2 cm 

3.6 cm 

Thermistor bolometer in conical light pipe 
1-2 ms 
10.5 sec 
20 - 80 cps 

Every second fringe of monochromatic 
source 5852.5 A 

325 samples sec" 1 

3408 

16 samples per interferogram (sampled 
4 times) 

3472 samples 

12 

41664 


Interferograms per orbit ~ 420 

Total number of bits per orbit 18 • 10 6 

Field of view 1.57 • 10' 2 ster(~8°) 

In-flight calibration after every 14 interferograms 

a. Outer space 

b. 300 °K blackbody 


35 



Monochromatic source 
Wavelength of reference source 
Detector for reference source 
Frequency of reference signal 


Neon lamps plus interference filter 
5852.5 A 
Silicon junction 
625 cps 
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